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The dynamics of expansion and large scale structure formation in the multicomponent Universe 
with dark energy modeled by the minimally coupled scalar field with generalized linear barotropic 
equation of state (EoS) are analyzed. It is shown that the past dynamics of expansion and future of 
the Universe - eternal accelerated expansion or turnaround and collapse - are completely defined by 
the current energy density of a scalar field and relation between its current and early EoS parameters. 
The clustering properties of such models of dark energy and their imprints in the power spectrum 
of matter density perturbations depend on the same relation and, additionally, on the "effective 
sound speed" of a scalar field, defined by its Lagrangian. It is concluded that such scalar fields 
with different values of these parameters are distinguishable in principle. This gives the possibility 
to constrain them by confronting the theoretical predictions with the corresponding observational 
data. For that we have used the 7-year WMAP data on CMB anisotropics, the Union2 dataset on 
Supernovae la and SDSS DR7 data on luminous red galaxies (LRG) space distribution. Using the 
Markov Chain Monte Carlo technique the marginalized posterior and mean likelihood distributions 
are computed for the scalar fields with two different Lagrangians: Klein-Gordon and Dirac-Born- 
Infeld ones. The properties of such scalar field models of dark energy with best fitting parameters 
and uncertainties of their determination are also analyzed in the paper. 
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I. INTRODUCTION 



The unknown nature of dark energy remains the most 
intriguing problem of physics and cosmology of the last 
decade. Theorists have proposed a dozen of models well 
matching the observational data (see e. g. reviews [l| and 
book [2]), but key tests for distinguishing of them still 
have to be developed. This can be achieved by the de- 
tailed study of their physical properties and confronting 
the theoretical predictions of each of them with obser- 
vations. Physically motivated and historically first after 
the cosmological constant dark energy candidate is the 
scalar field violating the strong energy condition. Here 
we analyze two such fields - classical and tachyonic - in 
order to find out how they match the current observa- 
tional data assuming that one of these fields dominates 
our Universe today. As in this paper we are interested 
mainly in the dark energy effects we restrict ourselves to 
models with zero curvature of 3-space. 

Field and fluid approaches are used for inclusion of 



the scalar field in the cosmological model, which de- 
scribes the dynamics of expansion of the Universe and 
formation of its large-scale structure. They are equiv- 
alent if the intrinsic entropy of dark energy is taken 
into account in the equations of perturbations. Here 
we use the combined approach in which the scalar field 
as dark energy is described by its dimensionless density, 
fide = 87rGp( de )/3iJoi the equation of state (EoS) pa- 
rameter w = P(de)/p(de) an d the effective sound speed 1 
c'j. = 5P(de)/v~P(de), determined by its Lagrangian. We 
assume also that a scalar field is minimally coupled to 
other components, i. e. it interacts with them only grav- 
itationally. Therefore these quantities completely define 
its model. In general case the last two values, w and 
e~, are arbitrary functions of time, which should be de- 
fined by additional assumptions about the properties of 
dark energy. The simplest case w = const and c s = 1 
(classical scalar field) is well studied. The recent deter- 
mination of value of w for this case by [H, based on the 
seven-year WMAP data (WMAP7), the distance mea- 



*Electronic address: novos@astro.franko.lviv.ua 

^Electronic address: olka@astro.franko.lviv.ua 

* Electronic address: apus@astro.franko.lviv.ua 

§ Electronic address: pelykh@iapmm.lviv.ua 



1 The terms "effective sound speed" and "adiabatic sound speed" 
of dark energy are used in the literature for designations of dark 
energy intrinsic values which formally correspond to thermody- 
namical ones. 
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surements from the baryon acoustic oscillations (BAO) 
in the space distribution of galaxies from SDSS DR7 
data @, Hubble constant measurements [|| and super- 
nova data pjj (WMAP7+BAO+SN) gives the value close 
to -1: w = -0.98 ± 0.053 (flat Universe, 68% CL). So, 
the cosmological model with such scalar field is rather 
similar to ACDM-model and has the same problems, e. 
g. fine tuning and cosmic coincidence. 

Varying in time EoS parameter w(a) of the so-called 
dynamical dark energy is the more complicated case 
which can in principle resolve those problems as well 
as promises the interesting properties of dark energy 
and its effect on the dynamics and structure formation 
of the Universe. The dynamical classical scalar fields 
called "quintessence" with different potentials were in- 
vestigated extensively during for over than a decade (for 
some early works see Q). Often the dependence of 
EoS parameter w on the scale factor a is assumed in 
ad hoc manner. The linear form of this dependence 
w(a) = wq + (1 — a)w a was proposed in [l2| and widely 
used. Here, wq and uu a denote the present values of w 
and its first derivative respectively. The determination of 
them by @ on the base of WMAP7+B AO+SN data gives 
w = -0.93 ± 0.12, w a = -0.41 ± 0.72. Such dark energy 
evolves from phantom field with w = —1.34 at the early 
epoch to the quintessential one at the current epoch with 
w = —0.93. The dark energy density increases from zero 
to ~ 1.25p° e at a « 0.83 and decreases asymptotically 
to zero after that. The previous determination of these 
parameters by [l3[ using similar but older datasets gave 
w = -1.09 ± 0.12, w a = 0.52 ± 0.46. It means that in 
this case the dark energy evolves from quintessential field 
with w = —0.57 in the early epoch to the phantom one 
at current epoch with w = —1.09. Its density decreases 
at early epoch, achieves the minimal value of ~ 0.83p^ e 
at the same scale factor a « 0.83 and grows later. Re- 
cently this form has been modified by [13j in order to 
bring the behavior of dynamical dark energy at early 
epoch closer to that of A-term. Their approximation 
has the additional third parameter a tran which, however, 
is weakly constrained by observations. The mentioned 
parametrizations of time dependence of EoS parameter 
allow the phantom divide crossing (w = —1) and extend 
the variety of properties of dark energy and its possi- 
ble physical interpretations. Other three, four and more 
parametric approximations of EoS have been analyzed 
[14| too. The additional degeneracies and uncertainties of 
parameters related to the early dark energy density and 
time variations of EoS parameter are inherent for them. 
And vice versa, the value of EoS w as well as of energy 
density f^ e related to the late epoch are determined well 
as a result of their main impact on the expansion his- 
tory of the Universe, horizon scale, distance to CMB last 
scattering surface and scale- independent growth factor of 
linear matter density perturbations. These values, how- 
ever, give no possibility to constrain essentially the types 
of cosmological scalar fields, or, in other words, the forms 
of their Lagrangians and potentials. 



The additional constraints on the type of cosmologi- 
cal scalar fields can be obtained by determination of the 
third mentioned above parameter - the effective sound 
speed cf, since it is related to the field Lagrangian as 

° 2 s = ml ( 2X m^ + is?)- lt governs the evolution of 
dark energy density perturbations which influence the 
evolution of matter density ones via the gravitational 
potential. Its impact on the linear power spectrum of 
matter density perturbations is essentially lower than the 
growth factor caused by the background dynamics, but it 
is scale-dependent and can be appreciable for some types 
of scalar fields. Some of them are studied carefully in 




In the previous papers [lj], Gjl we have analyzed the 
cosmological models with minimally coupled classical and 
tachyonic scalar fields with the EoS parameter varying 
in time as w(a) = woa 3 /(l + wo — woa 3 ) which was 
deduced from the condition of "zero adiabatic sound 
speed": c\ = P(d e )/P(de) = 0. It was shown that at 
the early epoch such fields can mimic the dark matter 
and diminish the fine tuning problem. The magnitude of 
dark energy perturbations of such fields increases before 
entering the acoustic horizon and decays in oscillating 
manner after that. It was shown also that they affect 
the magnitude of gravitational potential and cause the 
scale dependence of matter density perturbations which 
can be used for constraining the parameters of such dark 
energy models along with other cosmological ones. 

Here we analyze more general case - "constant adi- 
abatic sound speed" c\ = const, that is equivalent to 
assumption of the linear barotropic EoS of dark en- 
ergy P(de) — CiP(cie) + where C\ and Ci are con- 
stants. The important question related to such scalar 
fields is: how their introduction agrees with the obser- 
vational data on accelerated expansion of the Universe, 
its large scale structure and CMB anisotropy? To find 
the answer we perform analysis using the Markov Chain 
Monte Carlo (MCMC) method for multicomponent cos- 
mological model and the present data, namely seven-year 
WMAP data on CMB anisotropy SDSS DR7 and the 
SNIa photometric relations. 

The paper is organized as follows. In Section [H] we 
discuss the background dynamics of our models and 
parametrization of EoS parameter by the "adiabatic 
sound speed" . In Section IIIII the evolution equations for 
the scalar field perturbations in the synchronous gauge 
are presented. In Section IIVI the method of determina- 
tion of cosmological parameters and used observational 
datasets are described. In Section [V] we present the re- 
sults of determination of scalar field parameters and dis- 
cuss the properties of models with the best fitting ones. 
The conclusion can be found in Section IVI1 
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FIG. 1: Top panel: the dependences of EoS parameter on 
scale factor for barotropic dark energy with wo=-0.85 and 
different c 2 a (0, -0.2, -0.85, -0.99). Bottom panel: the depen- 
dences of dark energy density (in the units of critical one at 
the current epoch) on scale factor with EoS parameters pre- 
sented in the top panel. 



II. DARK ENERGY WITH BAROTROPIC EOS 
AND BACKGROUND EQUATIONS 

We assume that background Universe is spatially flat, 
homogeneous and isotropic with Friedmann-Robertson- 
Walker (FRW) metric of 4-space 



ds 2 



gijdx*dx j = a 2 (77) (dry 2 - 5 a pdx a dx p ), 



where r\ is the conformal time defined by relation 2 dt = 
a{7])drj and a{rf) is the scale factor, normalized to 1 at 
the current epoch rjQ. The Latin indices i, j, ... run from 
to 3 and the Greek ones are used for the spatial part 
of the metric: a, (3, ...=1, 2, 3. We also suppose that 
the Universe is filled with non-relativistic particles (cold 
dark matter and baryons), relativistic ones (thermal elec- 
tromagnetic radiation and massless neutrino) and mini- 



FIG. 2: The dynamics of expansion of the Universe with 
barotropic dark energy: H 2 (a) (top panel) and q(a) (bottom 
one) for the same models as in Fig. [T] 



mally coupled dark energy. The last one is considered 
as the scalar field with either Klein-Gordon (classical) or 
Dirac-Born-Infeld (tachyonic) Lagrangian 

L clas =X-U(cf>), L tach = -U(t)Vl - 2X, (1) 

where U(<p) and £/(£) are the field potentials defining the 
model of the scalar field, X — (j)-^ 1 /2 and X — £ ; i£ ;l /2 
are kinetic terms. We assume also that the background 
scalar fields are homogeneous (cj)(x,r]) — <fr(r)), C( x j^) = 
£(77)), so their energy density and pressure depend only 
on time: 



Pclas ' 
Ptach 



\J\-2X 



x - u {<!>), 



(2) 



The EoS parameters Wde 
following: 



P tach = -U(0\/l-2X. (3) 
= Pde/ Pde for these fields are 



Wclas = 



x-u 
x + u' 



Wtach =2X -1. 



(4) 



Here and below we put c = 1. 



Using the last relations the field variables and potentials 
can be presented in terms of densities and EoS parame- 
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ters as: 



4>{a) - 4> Q = ± 
U(a) = 



a da 1 ^ p de {a'){l + w{a')) 
a'H(a') 

p de (a) [1 - w(a)] 



for the classical Lagrangian and 



£(a)-e„ = ±y 



a'H(a') : 



(5) 
(6) 

(7) 
(8) 



for the tachyonic one. 

The dynamics of expansion of the Universe is com- 
pletely described by the Einstein equations 



where R, 



is the Ricci tensor and t\]\ T^ e 



(9) 



: the 



energy-momentum tensors of non-relativistic matter (m) , 
relativistic matter (r), and dark energy, (de), correspond- 
ingly. Assuming that the interaction between these com- 
ponents is only gravitational, each of them should satisfy 
the differential energy-momentum conservation law sep- 
arately: 



rpi (n) _ „ 

3 \i 



(10) 



Hereafter ";" denotes the covariant derivative with re- 
spect to the coordinate with given index in the space 
with metric gij and (n) stands for (m), (r) or (de). 

For the perfect fluid with density p/ n \ and pressure 
P( n )i related by the equation of state P(„) = W(„)P(„), it 
gives 



P(n) 



"3-P(n)(l + w (n)), 

a 



(11) 



here and below a dot denotes the derivative with re- 
spect to the conformal time, "'"= d/drj. For the non- 
relativistic matter W( m ) = and P( m ) = P(H) a ~ 3 ' l° r the 
relativistic one it>( r ) = 1/3 and p^ — p|°ja~ 4 . Here- 
after "0" denotes the current values. The EoS parameter 
W(de) and adiabatic sound speed c 2 ^ = P(de)l P(de) are 
related by the ordinary differential equation: 



Za- l (l+w){w-cl), 



(12) 



where a prime denotes the derivative with respect to the 
scale factor a. Here and below we omit index de for 



and 



'a (de) 



As it can be easily seen, the derivative 



of EoS parameter with respect to the scale factor will 
be negative for w < c 2 and positive for w > c 2 . In 
the first case the repulsive properties of scalar fields will 
be raising, in the second one - receding. In the general 
case c 2 can be an arbitrary function of time, but here 
we assume that it is constant: c 2 = const. In such case 



the temporal derivative of P(d e ) (v) 1S proportional to the 
temporal derivative of P(de)(v)- The integral form of this 
condition is the generalized linear barotropic equation of 
state 



P(de) = 



'aP(de) 



(13) 



where C is a constant. Cosmological scenarios for the 
Universe filled with the fluid with such EoS 3 have been 
analyzed in [2l[ . The solution of the differential equation 
(Tr2"l) for c\ — const is following: 



(1 + (£)(! + «*,) 



i + w a - (too - c 2) a 3(i+a 



1. 



(14) 



where the integration constant of (|12[) wo is chosen as 
the current value of w. One can easily find that (fT4")l 



gives (TT3|) with C 



(o) / 

P(de)( W 



;), where is cur- 
rent density of dark energy. Thus, we have two values 
Wo and c a defining the EoS parameter w at any red- 
shift z = aT 1 — \. The time dependences of barotropic 
EoS parameter for different values of c\ are shown in 
the top panel of Fig. [T] As it follows from (TT4")) . c\ 
corresponds to the EoS parameter at the beginning of 
expansion (u>,- m -t = w(0) — c 2 a , a = 0, z — oo). 

The differential equation (TTTI) with w from (TT?]) has the 
analytic solution too: 



(o) 

P(de) ~ P {de ) 



(1 + Wo )a- 3(1+c2j + c 2 a 



w 



1 



(15) 



The dependences of dark energy density on scale factor 
for different values of c\ are shown in the bottom panel 
of Fig. CD 

Using the obtained dependences of densities of each 
component on the scale factor one can deduce from the 
Einstein equations © the following equations for back- 
ground dynamics: 



H 

9. 



Ho^VLr/a^ + n m /a 3 + VL k /a 2 + n de f(a), (16) 

1 2fi r /a 4 + O ro /a 3 + (1 + 3w)il de f{a) 

2 n r /a 4 + n m /a 3 + n k /a 2 + n de f(a) 



(17) 



where /(a) = [(1 + w )a-^ 1+c ^ +c 2 a - w ]/(l + c 2 a ) and 
Clk = — K I (aH) 2 is the dimensionless curvature parame- 
ter. H = a/a 2 is the Hubble parameter (expansion rate) 
for any moment of time and q = — (ad /a 2 — l) is the 
acceleration parameter. They completely describe the 
dynamics of expansion of the homogeneous isotropic Uni- 
verse. 

To analyze the properties of scalar field models of dark 
energy with barotropic EoS and to calculate the observa- 
tional predictions the allowed by physics and mathemat- 
ics ranges of values for wq and c 2 should be determined. 



Often called in literature "wet dark fluid" . 
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FIG. 3: Evolution of potentials and kinetic terms of classical (left) and tachyonic (right) scalar fields with barotropic EoS for 
(? a > u>o (top panels), <? a — wo (middle panels) and c? a > wq (bottom panels). In the top panels = —0.2, in the middle ones 
Ca = —0.85 and in the bottom ones c„ = —0.99. In all panels wo = —0.85. 



The field variables 4> (f° r positive energy density of the 
scalar field) and £ are always real if w > —1. The poten- 
tial U (a) of classical scalar field is positive if w < 1 and 
the potential U(a) of tachyonic one is real if w < 0. As 
w(a) changes from c 2 to wq in the range of < a < 1, 
the allowed ranges of values for wq and c 2 are the same 
for both fields. Other upper constraint for classical scalar 
field can be obtained from the next speculations. In or- 
der to keep the shape of angular power spectrum of CMB 
temperature fluctuations close to the ACDM one match- 
ing well the observational data, the dimensionless density 
of dark energy Slde{a) — 8nGpd e (a)/3H 2 (a) should not 
exceed 0.01-0.02 at the last scattering surface. Its asymp- 
totic at small a is n de (a) w ^ e ) (l+w )a- 3 ( 1+c 'V(l+c? I ). 
The condition tt de (0.001) < 0.02 is satisfied for c\ < 0. 
Therefore, for both quantities wq and c 2 we accept the 
ranges of values — 1 < wq, c 2 < for which the scalar 
field variables and potentials have no singularities during 
the past history of the Universe, but there is no warranty 
that the singularities will not appear in future. Let us 
analyze that. 



The system of equations p^|) -([T7 1) describes 3 possible 
variants of future evolution of the Universe defined by 
the relationship between wq and (? a for given values of 
rest of the parameters f2 m , f2fc, fi r and Vide (Fig- D32]). 

1) w' < (c 2 a > w ): As it follows from (flip, in 
this case w decreases monotonically from c 2 at the early 
epoch to wq at current one up to -1 at the infinite 
time. The constant C in EoS (fT3")) is negative. The 
dark energy density and pressure tend asymptotically to 

P d 7 } = p£?(«£-«*)/(l + c2) and Pt } = -P^- There- 
fore, in this case the scalar field rolls slowly to the vacuum 
(see Fig. [1} and in far future the Universe will proceed 
into de Sitter stage of its expansion with q(°°) = — 1 and 
= V^e(c;i - w )/{c 2 a + 1)H . The scalar fields of 
such type have the following general properties (see rela- 
tions ((U)-© and top panels of Fig. [3]): their kinetic terms 
and potentials have always real positive values; for clas- 
sical field X < U and for tachyonic one X < 1/2; their 
potentials roll slowly to minima (U m i n , U m i n > 0), which 
correspond to the infinite values of the field variables and 
time; the kinetic terms of these fields tend asymptotically 
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FIG. 4: The w — dw/dlna phase plane for scalar fields with barotropic EoS as models of dynamical dark energy (solid lines). 
If dw/dlna < the fields evolve from right to left raising their repulsion properties, if dw/dlna > the fields evolve from left 
to right receding them. Thick dashed lines show the ranges occupied by the thawing and freezing scalar fields deduced by [12] 
from the analysis of simplest particle physics scalar field models of dynamical dark energy. In the left panel the phase plane 
evolution tracks of the scalar fields with barotropic EoS are shown in the range 0.5 <a<l(0<z<l) and in the same 
scale as in [22| for easy comparison. In the right panel the phase plane evolution tracks of the scalar fields with barotropic 
EoS correspond to the range 0.0001 < a < 10 (—0.9 < z < 10000). Thick solid black lines show the limits for such scalar field 
models: the upper line corresponds to c„ = — 1, the lower one to (? a = (superimposed with the lower limit for freezing scalar 
fields from |22|). The blue solid lines show the phase tracks of models considered in this section, the red solid lines show the 
phase tracks of the best fitting models discussed in section [V] The red point shows the best fitting w — const model. 



to - this means that </> and £ — > and the fields will 
"freeze" . 

2) w' = (c 2 = wo): It corresponds to the well- 
studied case w = const. In this case C = and we 
have usual barotropic EoS Pde = woPde, Pde — > when 
a — > oo. So, the future Universe will experience the 
power law expansion with a oc f 2 / 3 ( 1 +' t "o) anc j accelera- 
tion parameter q — > (1 + 3wo)/2. The scalar fields of 
such type have the following general properties (see rela- 
tions ((U)-© and middle panels of Fig. [3]): their kinetic 
terms and potentials have always real positive values; for 
classical field U /X = const > 1 and for tachyonic one 
X = (1 + wq)/2 = const; their potentials roll very slowly 
to minima, which correspond to the infinite values of the 
field variables and time; the potentials of these fields tend 
asymptotically to - it means that such fields "thaw" . 

3) w' > (c 2 < wq): The EoS parameter w 
increases monotonically from c 2 at the early epoch 
to wo at the current one and still continues to in- 
crease after that. It will reach in future at 

i 

a {w= o) = [cl(l + w )/ (cl - w )] 3 < 1 + c a) when poten- 
tial of tachyonic scalar field becomes imaginary (re- 
lation ([8]) and right bottom panel of Fig. [3]) and 

1 at a [w=l) = [i(l-cl)(l+w y(w -cl)]^^ 



when potential of the classical one becomes nega- 
tive (relation and left bottom panel of Fig. [3]). 
The scalar field densities at these a are positive: 

Pde(a( w=0 )) = P d °e(cl - wo)/c 2 a and p d e{a{w=x)) = 

Pd°e( c a ~ w o)/( c a ~ 1) correspondingly. The field will 
satisfy the strong energy condition p^ e + iPde > start- 
ing from a = [(1 + w )(l + 3c 2 a )/(2(c 2 a - w ))] 
and then accelerated expansion of the Universe will 
be changed by the decelerated one. The density of 

scalar field continues decreasing, reaches at a^ p=0 ) = 
i 

[(1 + wo)/(wo — Cq)] 3<1+c »' and then becomes negative. 
The EoS parameter at this moment will have discontinu- 
ity of the second kind (Fig. [TJ. Later, when p m + pd e 
reaches 0, the expansion of the Universe will be changed 
by the contraction since at this moment a = 0, a < 0, as 
it follows from equations fp~6|) and ([l~7l) , which have no so- 
lution for larger a. Such behavior can be corrected only 
slightly by the curvature parameter from the allowable 
range. Though the potential and field variable of tachy- 
onic scalar field become imaginary at a( w=0 ) < a < a( p=0 ) 
and a > a-( P =o) correspondingly (their energy density and 
pressure are real value always) we have no substantial ar- 
guments for ruling out such models since they have not 
singularities in past and are not excluded by the observa- 
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tional data. We have also no physical argument against 
the possibility of generation of such fields at the early 
epoch. Meanwhile, we will use the condition c 2 > wq 
too for analysis of the scalar field dark energy without 
singularities in future. The appearance of imaginary val- 
ues of either field variable or potential of the tachyonic 
scalar field suggests that this case needs to be analyzed 
using the multifield model of dark energy. That will be 
the matter of a separate paper. 



At the end of the analysis of general properties of the 
scalar field models of dark energy with c 2 = const let 
us examine their occupancy of the w — dw/dlna phase 
plane. From (|T^|) and the constraint — 1 < c 2 < fol- 
lows that the scalar field models of dark energy with 
c 2 = const occupy the w — dw / din a region limited by the 
lines dw/dlna = 3(f + w) 2 and dw/dlna = 3w(l + w) 
(Fig. 0}. The last one coincides with the lower limit 
for freezing scalar field models of dark energy deduced 
by (22| from the analysis of the simplest particle-physics 
models of cosmological scalar fields. Below it the scalar 
fields have too large density at the early epoch that con- 
tradicts the data on CMB anisotropy. Above the upper 
limit there is a range of fields which started as phantom 
ones, which is excluded for fields with classical and tachy- 
onic Lagrangians considered here. The scalar fields which 
are in the phase plane between the lines dw/dlna = 
and dw/dlna = 3w(l + w) evolve from right to left in 
Fig. [4] and their repulsion properties are raising with 
time. They are unlimited in time and w for them tends 
asymptotically to -f. The scalar fields which are in 
the phase plane between the lines dw/dlna = and 
dw/dlna = 3(1 + w) 2 evolve from left to right in the 
Fig. 0] and their repulsion properties are receding with 
time (dw/dlna > 0, w increases) to change the accel- 
erated expansion by decelerated one and even collapse. 
They can start in the range below the lower limit for 
thawing scalar fields, then enter the range of thawing 
scalar fields limited by [13, cross it and go out of up- 
per limit dw/dlna = 3(1 + w) when w > 0. So, the 
scalar fields with dw/dlna > (c 2 < wq) can only par- 
tially be called thawing. We propose to call them "scalar 
fields receding repulsion", reflecting their main proper- 
ties. Symmetrically, the scalar fields with dw/dlna < 
(c 2 > wq), occupying the same range as freezing scalar 
fields from [13 , can be called "scalar fields raising repul- 
sion" . Establishing of how well they can fit the current 
observational data is the goal of the paper. 



Since the classical scalar fields are completely indistin- 
guishable from the tachyonic ones when considering only 
the background dynamics, the linear theory of perturba- 
tions in the multicomponent Universe must be included 
into analysis. 



III. EVOLUTION OF LINEAR 
PERTURBATIONS IN THE SYNCHRONOUS 
GAUGE 

Now we consider the spatially flat (fife = 0) back- 
ground with scalar perturbations in the synchronous 
gauge with the line element 

ds 2 = g 2J dx l dx j = a 2 (n)(dn 2 - (6 af) + h afi )dx a dx [i ) (18) 

The scalar perturbations of metric h a p can be decom- 
posed into the trace h = and traceless h a p compo- 
nents as h a/ 3 = hda/3/3 + h a/ 3. The perturbations are 
supposed to be small (h <C 1), henceforth all following 
equations are linearized with respect to the perturbed 
variables. In the multicomponent fluid each component 
moves with a small peculiar velocity V a — dx a /dn, de- 
fined by its intrinsic properties (density, pressure, en- 
tropy etc.) and h. At the linear stage of evolution of 
perturbations the cold dark matter (CDM) component is 
a pressurelcss perfect fluid interacting with other compo- 
nents only via gravity. Therefore, the synchronous coor- 
dinates are usually defined as comoving to the particles of 
CDM. The energy density /?(„) and pressure P(„) of each 
component are perturbed as p(„) (a; Q , 77) = P(„)(?/)(l + 
$(n)(x a ,r])) and P {n) (x a ,rj) = P (n) (n) + SP (n) {x a ,n). 
They are measured by a comoving observer being at rest 
relative to the fluid at the instant of measurements, so 
their perturbations are related as <5P(„) = c^ n ^p( n )^(n), 
where c 2 ,s is the comoving effective sound speed of com- 
ponent (n). This sound speed equals to zero for dust 
medium (e. g. CDM or baryonic matter after recombi- 
nation at scales above Jeans scale) and 1 /3 for relativistic 
components. For scalar fields the entropy perturbations 
are inherent and cause in addition to the adiabatic pres- 
sure perturbations, which follow from the variation of 
(fT3f , the non-adiabatic pressure ones SP^^ , so the to- 
tal perturbation is their sum [23j : 

SP(de) = c 2 a 5 P{de) +5Pl™ d) - 

In the rest frame of dark energy (V^e) = 8(f) = S£ = 0) it 
can be presented as 

SP(de) = c l(de)P{de)5(de)i 

where the effective (rest-frame) sound speed c 2 ^ for the 
scalar field with given Lagrangian can be calculated as 

2 _ PjX = L, x 
c s(de)- p x 2XL, XX + L.X 

It equals 1 for the classical scalar field and — w for the 
tachyonic one. 

For pure barotropic fluid 4 SP^^ = and c 2 — c 2 



4 The conditions on the Lagrangian for vanishing non-adiabatic 
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FIG. 5: Evolution of linear density perturbations of cold dark matter (dashed line), baryons (dotted) and scalar field (solid) 
with = wo = —0.85, (? a = 0, <? a — —0.2, e„ — —0.99 (from top to bottom). Classical scalar field - left column, tachyonic one 
- right. The wave number of perturbations is k — 0.1 Mpc -1 . 



- that distincts it from the scalar fields with barotropic 
equation of state. 

In the linear perturbation theory it is convenient 
to perform the Fourier transformation of all spatially- 
dependent variables and use the equations for corre- 
sponding Fourier amplitudes. The differential energy- 
momentum conservation law 5Tj .j = for perturba- 
tions in the space with metric (|18|) gives the equations 
for evolution of Fourier amplitudes of density and ve- 
locity perturbations of dark energy in the synchronous 
gauge: 

$(de) + 3(c^ - w)aH5 {de) + (1 + w)^ 



+(l + w) 



k + 9a 2 H 2 



V, 



(de) 



(19) 



V (de) + aH(l - 3c 2 s )V (de) - Y^ 6 iM) = °' ( 20 ) 
These equations in the conformal-Newtonian gauge are 



pressure perturbation in the case of generic scalar field can be 
found in [Hi. 



presented in [15[ , in the gauge-invariant variables in [19j . 
In these papers it was shown also that the perturba- 
tions of scalar fields with different c 2 evolve differently. 
For the case of classical and tachyonic scalar fields with 
c 2 = the evolution was extensively studied for the two- 
component (dark energy plus dark matter) model in [T|| . 
It was shown also that the transfer function of the cold 
dark matter perturbations is different for the same set of 
cosmological parameters but different Lagrangians (and 
thus c 2 s ) of dark energy. 

The equations for the rest of components^ (non- 

rclativistic and relativistic) are the same as in [25j, so 
we skip this part here. For calculation of the evolution 
of perturbations in all components and power spectra 
of matter density perturbations and cosmic microwave 
anisotropy we have used the publicly available code 
C AMB [27|, l28j , modified to include the presented here 
expressions for H(a), /3(de)(«) and the evolution equa- 
tions for dark energy perturbations. The evolution of 
density perturbations (A;=0.1 Mpc -1 ) of cold dark mat- 
ter, baryons and scalar fields with different c 2 (0, -0.2, 
-0.85, -0.99) is shown in Fig. O Other cosmological pa- 
rameters have the next values: Vt de = 0.722, wq = —0.85, 
fl cdm h 2 = 0.111, fl b h 2 = 0.0227, Hq = 66.1 km/s/Mpc, 
n s = 0.975, r re i — 0.085. The adiabatic initial conditions 
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FIG. 6: Top panels: The ratios of matter density power spectra P(k) in the models with classical (left) and tachyonic (right) 
scalar fields to the matter density power spectra in the ACDM model. Bottom panels: the ratios of CMB angular power spectra 
Ci in the models with classical (left) and tachyonic (right) scalar fields to the corresponding ones in the ACDM model. In 
each panel the solid line stands for the scalar field with c 2 a — wo — —0.85, dotted one for c„ = 0, dashed one for cl = -0.2, 
dash-dotted for c„ = —0.99. 



from 25j for all components except dark energy are used 
here and below. The initial conditions for dark energy 
perturbations are obtained from the asymptotic solutions 
of p^|) - (pU)l for fc?7 <C 1 (the super horizon perturbations) 
in the early radiation-dominated epoch and are following: 



t t init 

V (de) ~ 



(4-3c 2 )(l+ W ) 



8 + 6c 2 - 12w + 9c 2 s (w~cl) 



(? s kr]i, 



8 + 6c 2 s - 12w + 9c 2 s (w ~ c 2 a ) 



h, (21) 



h. (22) 



The character of evolution of scalar field density pertur- 
bations depends on the temporal behavior of EoS pa- 
rameter, the difference of adiabatic and effective sound 
speeds (the value of intrinsic entropy of scalar field is pro- 
portional to it) and the ratio of the scale of perturbation 
to particle horizon (see equations (fll?)) -([2T) |) ). We can see 
that the magnitude of density perturbations of such fields 
increases before entering the acoustic horizon, decreases 
and changes the sign during entering one and decays with 
oscillations later. So, at the current epoch the magnitude 



of dark energy linear density perturbations is essentially 
lower than matter one at all subhorizon scales. 

Evolution of matter density perturbations is affected 
by the dynamics of background expansion through the 
scale-independent expansion rate and by scalar field den- 
sity perturbations through the perturbations of metric 
which depend on scale. For visualization of such to- 
tal influence of dark energy on matter density evolution 
we present in Fig. [5] the ratios of matter density power 
spectrum P(k) and CMB angular power spectrum Ci in 
the models with scalar field to the corresponding ones 
in the ACDM model with the same cosmological param- 
eters. The dependence of both ratios on scale, c 2 and 
type of scalar field inspires hope that usage of the accu- 
rate enough dataset will give the possibility to constrain 
the scalar field models of dark energy considered here. 
For that we use here the available currently observational 
data. 

At the end of this section let us discuss the differ- 
ence between scalar field models of dark energy with 
barotropic equation of state and pure barotropic fluid 
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FIG. 7: Pure barotropic fluid models of dark energy in the multicomponent Universe. Left: the evolution of linear density 
perturbations (k — 0.1 Mpc -1 ) of cold dark matter (dashed line), baryons (dotted) and barotropic fluid (solid) with <? a — (? a = 
0.0, <? a = (? a — 0.1 and <? a — (? a — —0.1 (from top to bottom). Right: the ratios of matter density power spectrum P(k) (top 
panel) and CMB angular power spectrum Ci in the model with barotropic fluid to the corresponding ones in the ACDM model. 
In each panel the solid line corresponds to the pure barotropic fluid with c? s = (? a = 0.0, dashed one - with (? s — c? a = —0.1. 



models of dark energy. For pure barotropic fluid (PBF) 
c 2 = c 2 . The analysis of equation (49) from |l9| shows 
that such dark energy is strongly gravitationally instable 
at subhorizon scales of density perturbations (krj <C 1) for 
c 2 < 0. So, such fluid model of dark energy will not con- 
tradict the observational data when c 2 — c 2 > 0. In the 
left panel of Fig. [7] the evolution of linear density pertur- 
bations of cold dark matter, baryons and pure barotropic 
fluid with c 2 s = c\ = 0, c 2 s =c 2 a = 0.1 and c 2 = c 2 a = -0.1 
is shown for perturbations with the wave number k = 0.1 
Mpc -1 . The curves for last case (c 2 = —0.1) illustrate 
the strong instabilities in perturbation growth, so that 
the models with c 2 < should be avoided. In order 
to check the viability of PBF models of dark energy we 
have computed the ratios of matter density power spec- 
trum P(k) and CMB angular power spectrum C; in such 
models to the corresponding ones in the ACDM model. 
In the right panel of Fig. [7]one can see, that the influence 
of PBF dark energy with c 2 = c 2 = on the evolution of 
matter density perturbations differs from one of classical 
and tachyonic scalar fields with c 2 = (see Fig. [5]). The 
large difference between the power spectra of PBF+CDM 
with c 2 — 0.1 and ACDM models means that the upper 
observational constraint on c 2 should be close to zero. 



In the context of this discussion and results, presented 
in this section, we would like to explain the title of 
the analyzed here dark energy model - scalar field with 
barotropic equation of state. It was shown above that 
the scalar field is not exactly barotropic because of the 
non-adiabaticity, which appears at the level of linear per- 
turbations and is cx 5(de) ■ Since 5(d e ) *C 1 now and in the 
past at all scales (at large scales due to the initial power 
spectrum cx k n ° , at small scales due to the decaying of 
subhorizon perturbations of the scalar field, shown in Fig. 
[5]) , in the main order the dark energy equation of state 
is the generalized linear barotropic one. This property 
defines the scalar field dynamics as well as the dynamics 
of expansion of the Universe and its future. That is why 
we call such scalar field the "scalar field with barotropic 
equation of state" . 



IV. METHOD AND DATA 

To determine the best fitting values and confiden- 
tial ranges of the scalar field parameters together with 
other cosmological ones we perform the Markov Chain 
Monte Carlo (MCMC) analysis for the set of current ob- 
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FIG. 8: One-dimensional marginalized posteriors (blue lines) and mean likelihoods (black ones, averaged in each of 30 bins) 
for Q, de , w and <? a (from left to right) for the combined datasets WMAP7+SDSS LRG DR7 (upper row) and WMAP7+SDSS 
LRG DR7+SN+HST+BBN (bottom row). Solid lines correspond to the classical scalar field, dashed lines - to the tachyonic 
one. 



servational data, which include the power spectra from 
WMAP7 0] and SDSS DR7 0, the Hubble constant 
measurements Q, the light curves of SNIa [ll[ and Big 
Bang nucleosynthesis (BBN) prior (26j . 

We use the publicly available package CosmoMC 0, 
IH, which includes the code CAMB Hi [H| for calcu- 
lation of model predictions for sampled sets of 8 cos- 
mological parameters (for recent usage of this package 
see [3l|) : the present value of the dark energy EoS Wo, 
the physical energy density of cold dark matter il c h 2 , 
the physical energy density of baryons fl^h 2 , the Hubble 
constant 6 Hq, the scalar spectral index of the primor- 
dial perturbations n s , the amplitude of primordial per- 
turbations A s , the reionization optical depth r and the 
Sunyaev-Zel'dovich amplitude A$z- The run of original 
CosmoMC with dataset listed above gives the best fit- 
ting values and confidential ranges for these parameters 
which arc presented in the first and second columns of 
Table |H] 

The original CAMB has been modified as described 
above, so the CosmoMC has been run with proposed here 
parametrization of dark energy EoS parameter (|14[) . The 
extra parameter c 2 can be either fixed or left free for 
determination of its best fitting value. Thus, we have 



5 The curvature parameter was fixed at Qk = 0, so Qde = 1 — 

fifr ^cdm ■ 

6 Instead of Hq the CosmoMC code varies the parameter 9 which 
is the ratio of the sound horizon to the angular diameter distance 
multiplied by 100. 



9 free parameters to be explored with MCMC. By now, 
we have excluded the phantom models (wde < — 1) and 
assumed flat priors for wq and c 2 in the range —1 < 
W , c 2 a < 0. 



V. RESULTS AND DISCUSSION 

We have performed MCMC runs for 9 spatially flat 
cosmological models with the dataset described above. 
All models consist of 5 components: scalar field as dark 
energy, dark matter, baryons, thermal background ra- 
diation and massless active neutrinos. The models dif- 
fer by the type of scalar field (classical - CSF or tachy- 
onic - TSF) and additional conditions for values of adi- 
abatic sound speed (c 2 = Wq, c 2 = 0, c 2 > wq and 
c 2 free). We assume that dark matter is cold (CDM) 
and since the density of dark components are dominating 
now and defining the dynamical and clustering proper- 
ties of our Universe the models are called CSF+CDM or 
TSF+CDM with corresponding additional conditions for 
c 2 which determine the type of dark energy. Each run 
has 8 chains and the number of samples in each chain is 
~ 200 000. The main results are presented in the Tables 
U and [Hi where the parameter values of best fit samples 
and lcr confidential ranges are shown for the main cos- 
mological parameters. In the last row of the tables the 
quantity — logL = % 2 /2 of best fit sample is presented 
for each model. 

The most interesting models with free c 2 
we have tested using 2 sets of observational 
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TABLE I: The best fitting values and la confidential ranges of 
cosmological parameters in the CSF+CDM and TSF+CDM 
models, determined on the base of WMAP7 and SDSS LRG 
DR7 data. The current Hubble parameter Ho is in units 
kms~ 1 Mpc _1 , the age of the Universe to is given in Giga 
years. 



Parameters 


CSF+CDM 


TSF+CDM 


Wo 

cl 

WQcdmh 2 

Ho 

n s 

log(10 10 A s ) 
to 


o.69t° ;°o 6 8 

-0.801^ 

-0 Q1+ - 91 

9 9 7 +0.16 
-0.15 

i m +0 - 14 

65.4±™ 

o-97±g:SS 

3.08±g;Sg 
10.8±|;1 
n 7^+°- 13 

u -'°-0.21 

13.8t°; 6 3 


0.68ir o 6 8 
-0 20+ 19 

9 97+O.I7 
-0.16 

., -,9+0.12 
■ L " LZ -0.1S 

65.01^ 
0.97«;° 4 4 
3.08t° : °» 
10.4±|;f 
0.741^ 

u.o±°i 


-logL 


3760.26 


3760.39 



data: WMAP7+SDSS LRG DR7 (Table HJ) and 
WMAP7+SDSS LRG DR7+SN+HST+BBN (Table El 
last 2 columns). The normalized mean likelihood func- 
tions L{xj\9j) and marginalized posterior distributions 
L(9i]Xj) for fide, wo and (? a are presented for both 
datasets in Fig. HI The first of them, L(xj;9i), is the 
probability distribution of the observational data Xj for 
given model parameters 9i. The second one, L(xj]0i), 
is the probability distribution of the parameters 9i for 
given observed data Xj. 

We can see that for both fields the current density 
tide is constrained well by both dataset: the likelihoods 
and posteriors are close. The current density Qde and 
EoS parameters wq are determined more accurately by 
the complete dataset: fl de = 0.72^^;^, w = — 0.93±^of 
for the classical scalar field and fide — 0-72to'o5' w o = 
— 0.971q o3 for the tachyonic one. The adiabatic sound 
speed c„, which is the EoS parameter at early epoch, 
is essentially unconstrained: its la confidence range is 
wide and coincides practically with the prior range [-1, 



0]. The mean likelihood and posterior are different, for 
both fields the likelihood is bimodal. First peak is close 
to -1, another one to 0. For classical scalar field the bi- 
modal character of distribution is more appreciable. It 
means that the used datasets are not appropriate for the 
estimation of this parameter. On the other hand it means 
that for any fixed value of c 2 . from the range [-1, 0] there 
exists a set of the rest best fitting parameters for which 
the model predictions match the used observational data 
well. The results presented in the columns 1-7 of Table 
[II] illustrate that. In order to remove the uncertainties in 
determination of c 2 , the other data and/or other statis- 
tical methods of analysis must be used, but this will be 
the matter of a separate paper. 

Analogical likelihood and posterior distributions for 
the rest of parameters have the shape of the Gaussians 
with half width equal approximately to the values shown 
in the upper/lower indices of each value in the tables. We 
can see also that the best fitting values of all parameters 
of models presented in the Table |TT] in columns 1-7 lie 
in the la marginalized confidential limits of the models 
with free c 2 . The differences of x 2 values are statistically 
insignificant. It means that observational datasets used 
here give no possibility to distinguish these 9 models at 
high statistical level, none of which can be ruled out by 
these data. 

In Fig. [9] we compare the power spectra of CMB tem- 
perature fluctuations and matter density ones for cos- 
mological models with the best fitting parameters from 
Table HI with WMAP7 and SDSS DR7 LRG data. They 
demonstrate perfect agreement of all models with obser- 
vations. The power spectra of best fitting ACDM model 
(its parameters are presented in the column 1 of Table HI)) 
overlap with the spectra of CSF+CDM and TSF+CDM 
models. 

Let us discuss the best fitting CSF+CDM and 
TSF+CDM models presented in the Table HH For com- 
parison we present the ACDM model with the best fitting 
parameters determined for the same dataset by original 
CosmoMC. 

1. Scalar field models with c 2 = (columns 4 and 
5 in Table HI]) . Such fields have partially the properties 
of ACDM model. Their pressure is constant during the 
whole history of the Universe: Pde — w aPd°e ~ cons t- 
But their density pde — Pdeli^ + w o) a ~ 3 — wq] at the 
early epoch changes similarly to one of the dust mat- 
ter (oc a -3 ) and asymptotically approaches the constant 

value — lOopSe m ^ ar future (a 3> 1). The EoS parame- 
ter w(a) = wga 3 /(l + wq — woa 3 ) decreases from m at 
the early epoch (a w 0) to -1 in future. Now it is very 
close to its asymptotic value: wo ps —0.99. The expan- 
sion rate ([TrJ]) decreases from H(a) = y/TE^Hoa~ 2 in early 
epoch to the constant value V. = \Z—w QdeHo ~ 0.85H 
in future. The acceleration parameter varies slowly from 
+1 in the early Universe to = —0.57 for classical field 
and go = —0.58 for tachyonic one at the current epoch 
up to -1 in far future. So, such scalar fields in the early 
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TABLE II: The best fitting values and la confidential ranges of cosmological parameters in the ACDM, CSF+CDM and 
TSF+CDM models determined on the base of complete data set listed in the section [IV] The current Hubble parameter Ho is 
in units kms -1 Mpc -1 , the age of the Universe to is given in Giga years. 



Parameters 


ACDM 
1 


CSF+CDM 
w=const 
2 


TSF+CDM 
iu=const 
3 


CSF+CDM 
cl = 
4 


TSF+CDM 

cl = 
5 


CSF+CDM 

cl > Wo 
6 


TSF+CDM 
cl > Wo 
7 


CSF+CDM 
cl free 
8 


TSF+CDM 
cl free 
9 


fide 


n 79+ 04 

u - ' z -0.05 


72+ ' 04 

u - ' z -0.04 


72+ ' 04 

u - ' z -0.05 


72+ ' 04 

u - ' z -0.05 


72+ 04 

u - ' z -0.05 


71 + 04 

u -' i -0.06 


u - ' 1 -0.04 


72+ 04 

u - ' z -0.06 


72+ - 04 

u - < z -0.05 


W 


-1 


-1 00+ 017 
r.uu_ 00 


_n qq+0-16 


_n qq+0-03 
u - yy -o.oi 


r.uu_ 00 


-0.99l°; « 


-i-oolS:S 


-0.93±°; 3 7 


-0.971°;^ 


cl 


-1 


1 00+°' 17 
-x.uu_ 00 


u - 3M -0 .01 








-0.051™! 


-o.o6in 6 4 


-o.97l° ;^ 


-0.99l° ; 9 o? 


100fi b /t 2 


9 07 +0.13 
-0.14 


9 or+0.16 
Z - ZJ -0.12 


2 26+ 15 

Z - ZD -0.15 


2.26to';i 4 


2 24+ - 17 

Z - Z4 -0.13 


9 9 r+0.17 
z - zd -0.12 


9 9 r+0.17 
Z - ZJ -0.13 


9 9 r+0.17 
z - zd -0.13 


2 29+ 14 

z - z »_0.16 


Wflcdmh 2 


l-14±oTs 


i i p+0.09 

1 - io -o.n 


1.141°;?? 


1.14lg;» 


i-i4^;i°3 


-i -lo+O.lO 
J— LO -0.13 


I ir+0.08 
i - 1O -0.14 


i 11+0.13 
1 - li -0.12 


i 1O+0.10 
J - xo — 0.13 


H 


7o.o±l : g 


69.4t« 


69.7I 4 4 ° 


69.2±|;| 


69.5l 4 ;° 


68.61H 


69.2l 4 ; 2 7 


69.21 4 ; 2 


69.71 3 ; 7 


n s 


qv+0-03 

<->•»( -0.03 




n 07+0.03 

-0.04 


0.98l :°oS 


Q7+0-03 
u - y ' -0.03 


-0.04 


Q7+0-04 
-0.03 


Q7+0-05 
u - y ' -0.03 


n Q7+0-04 
u - a ' -0.04 


log(10 10 A 3 ) 


3.09i° ;°o 9 9 


o nn+0 09 
O,u3 -0.09 


3.091°:°? 


3-09±° :io 


3.09l ;°09 


O QQ + 0.10 

°- ua -0.09 


3.09l°;J° 


3-071°;^ 


3.09l°; ° 




10.5l^ 3 


10.61^4 


10.5± 3 ; 2 


10.71 3 ; 7 


10.51 3 ; 2 


10 - 8 ±3.7 


10.31 3 ; 4 


io.3l 3 ; 9 3 


10.41 3 ; 4 


0"8 


n oq+o.06 


0.83lg;» 


0.83lH 6 


0.78t°;i° 


0.82«;? 7 


0.78l°;i° 8 


0.83t°;°j 


0.79l°;i» 


n 01+O.O8 


to 


o 7+0.3 
iJ - ' -0.3 


13.8±g-J 


i o 7 +0.3 
' -0.3 


io 7 +0.4 
10 - ' -0.3 


13.8l°; 3 


l3.8l° ;l 


13.8lg; 4 


13.8l°o1 


13 7 +0.4 
10 - ' -0.3 


-logL 


4027.60 


4027.63 


4027.62 


4027.24 


4027.75 


4027.51 


4027.79 


4027.35 


4027.41 



Universe mimic the dust matter and will mimic vacuum 
energy density in future. The final stage of evolution of 
the Universe with classical or tachyonic scalar field with 
cl = will be de Sitter one: a oc expHt forever. These 
fields jointly with cold dark matter can be considered 
as some type of the generalized dark component which 
decays and has properties of the two-component model, 
dark matter plus dark energy. The evolution of density 
perturbations of matter and both scalar fields is shown 
for the wave number k = 0.1 Mpc -1 in Fig. [5] (for other 
scales in the two-component model see The den- 

sity perturbations of classical and tachyonic scalar fields 
evolve differently because of difference of their effective 
sound speeds (c 2 = 1 for classical field and c 2 = —w(a) 
for the tachyonic one). Their impact on the matter den- 
sity perturbations in the models with the same parame- 
ters is distinguishable too (see Fig. [fj]). However, small 
variation of the main cosmological parameters (columns 
3 and 4 in the Table |ll| reduce this distinction (Fig. [9]) 
at the cost of a bit larger \ 2 for the tachyonic field. At 
the decoupling (z^ec ~ 10 3 ) and at the earlier epochs the 



ratio of scalar field density to matter one was constant 
{pde/Pm ~ (l + wo)ride/^m) and for parameters from the 
table it equals 0.015. So, scalar fields with such pa- 
rameters do not change practically the predictions of the 
concordance ACDM model for the Big Bang nucleosyn- 
thesis, however, they reduce essentially the fine tuning 
problem: at the end of phase transitions (i ~ 10~ 10 s) 
their energy densities were about 10 -12 of radiation en- 
ergy density against 10~ 54 for ACDM-model. 

2. Scalar field models with c 2 > wq (columns 6 and 
7 in Table [TT| • The condition c 2 > wq changes slightly 
the distributions shown in the bottom row of Fig. [8] but 
causes finding of the best fitting value of c 2 near 0, there- 
fore the properties of such dark energy models are sim- 
ilar to ones with fixed c 2 = 0. The future of the Uni- 
verse with such fields is de Sitter expansion a oc expUt, 

~ for 



«*>)/(! 



)H 



0.84ff, 



where H = yJVl de (c 2 a 
both fields. 

3. Scalar field models with w — const (columns 2 and 
3 in Table [TTJ) . The density of scalar fields with the best 
fitting parameters decreases very slowly {pde a 



i) 



-0.0003 



14 





FIG. 9: The power spectra of CMB temperature fluctuations (top) and matter density ones (bottom) for cosmological models 
with classical (left) and tachyonic (right) scalar fields and best fitting parameters from TablellTl The corresponding observational 
data from WMAP7 and SDSS LRG DR7 are presented. 



for classical field and oc a -0 03 for tachyonic one) during 
the whole history of the Universe, the pressure follows 
it rigorously: Pde ~ ~Pde- The ratio of scalar field den- 
sity to matter one was pde/ Pm ~ at the decoupling 
and lower at the earlier epoch. So, such fields also do 
not change the expansion dynamics of the early Universe 
and predictions of the concordance ACDM model. In the 
very early Universe, at t ~ 10 -10 s, the energy densities 
of classical and tachyonic scalar fields with the best fit- 
ting parameters were about 10~ 54 of radiation one, so, 
the models with such fields have the same fine tuning- 
problem. The acceleration parameters at current epoch 
are the same for both fields: qo = —0.57. Therefore, the 
past and future dynamics of expansion of the Universe 
with such fields is similar to the ACDM model one. The 
best fitting parameters of CSF+CDM and TSF+CDM 
models with w — const are close, since the evolution of 
their density perturbations is similar for values of w rel- 
atively close to -1 (see Fig. [5]and[n]). So, such fields with 
different Lagrangians are practically indistinguishable by 
cosmological observations. 

4. Scalar field models with — 1 < c\ < wq (columns 
8 and 9 in Table HI)) . The best fitting values of c\ and 
wq are -0.97 and -0.93 for classical scalar field, -0.99 and 
-0.97 for tachyonic one correspondingly. The EoS pa- 



rameter w increases very slowly from the first value to 
the second one during the age of the Universe. In the 
early epoch the density of classical scalar field changed 
oc a -0-09 , the density of tachyonic one oc a -0 04 . Their 
energy densities were only ~ 10~ 51 of radiation one at the 
baryogenesis epoch (t ~ 10~ 10 ) and ~ 10~ 31 at the nu- 
cleosynthesis epoch. At the decoupling they were ~ 10~ 8 
of matter density. For the ACDM model these values are 
~ 10~ 54 , ~ 10~ 32 and ~ 10~ 8 respectively. So, in this 
case - like the previous one - the classical and tachy- 
onic scalar fields do not change the expansion dynamics 
of the Universe in its past history and do not affect the 
local physical processes such as cosmological recombina- 
tion, nucleosynthesis, baryosynthesis etc. comparing to 
the concordance ACDM model. The current acceleration 
parameters are qa = —0.51 in the model with classical 
scalar field and qo = —0.55 in the model with tachyonic 
one. But predicted by these models future dynamics of 
the Universe differs essentially from the one predicted by 
previous models. The general properties were discussed 
in section [Til here we particularize them for the best fit- 
ting parameters of the model with classical scalar field. 
So, in this model the EoS parameter will continue to in- 
crease in future. It will reach -1/3 at ai q —Q\ ~ 530, when 
accelerated expansion of the Universe is changed by the 



15 



decelerated one (the age of the Universe will be pa 195 
Gyrs). Thereafter w will reach at ar w —Q\ ps 678 when 
the age of the Universe is pa 208 Gyrs. If the scalar field 
is tachyonic then at this moment its potential will be- 
come imaginary. The potential of classical scalar field 
will become negative when w reaches 1 at dw—i pa 804, 
the age of the Universe will be pa 226 Gyrs then. The 
next important moment in the evolution of the scalar 
field with such parameters will occur when the Universe is 
pa 268 Gyrs old (a = 950): it will violate the weak energy 
condition (pd e > 0) and its energy density will become 
negative after that. The pressure at this time will have 
constant positive value (Pde ~ {wo ~ c 2 ) pa 0.03/?^*), 
EoS parameter will have a discontinuity of the second 
kind (like in Fig. [T]), as well as the potentials of both 
fields and kinetic terms of tachyonic scalar field (like in 
Fig. [3]). Some time after a( Pdc=0 ) the total energy den- 
sity remains positive (p m + pde > 0), the strong energy 
condition p m + Pde + 3Pde > is still satisfied. The decel- 
erated expansion of the Universe will stop finally when 
the scale factor reaches the maximal value a max which 
is slightly larger than ai pde —Q\. At this moment a = 
and & < 0. The values of matter and dark energy den- 
sities at the turnaround are p m — —pde ~ 3 • 10~ w p^. 
The age of the Universe will be pa 268 Gyrs then. After 
that the collapse begins. So, the Universe, filled with the 
scalar field with such parameters, is limited in time by 
the age of pa 540 Gyrs and finishes its existence in the 
Big Crunch singularity. Of course, the possibility of exis- 
tence of scalar fields with such unusual properties needs 
yet the comprehensive analysis. 

Therefore, all models presented in Table [TT1 match well 
the set of current observational data used here (Fig. [SJ. 
Using different subsets (for example, excluding SDSS 
LRG DR7, SN data or both) or substituting the last re- 
lease of SDSS LRG DR7 by one of the previous ones, 
e. g. by DR4, changes somewhat the best fitting val- 
ues, so the preferable model has c a > wq. The same 
relation between these parameters holds when the SDSS 
date are excluded as well as WMAP7 and Union2 SN 
data are substituted by the WMAP5 [H and Union Q 
datasets correspondingly (see also [2(|). It means that 
the different datasets, obtained by different techniques 
and corresponding to different redshifts, are not in com- 
plete agreement with each other yet. 

VI. CONCLUSION 

We have analyzed the minimally coupled cosmological 
scalar fields with the time-variable EoS parameter w(a) 
and constant adiabatic sound speed c a , which means that 
the temporal derivative of pressure follows rigorously the 
temporal derivative of energy density: Pde oc Pde- In this 
case the equation of state has two natural parameters, wq 
and c a , which are the EoS parameter at the present epoch 
and at the early stages correspondingly. The evolution of 
w(a) from c 2 at the beginning to wq now is monotonous, 



however the future character of it is determined by the 
relation between the values of wq and c 2 . So, in the case 
c a > wq the derivative of w with respect to the scale 
factor a is negative, w' < 0, and fields roll down to the 
vacuum energy at the infinite time like one with w = — 1 . 
In the opposite case c 2 < wq, w' > and fields lose their 
repulsion properties. The accelerated expansion of the 
Universe will change in future to the decelerated one, the 
weak energy condition will be violated and the Universe 
will eventually collapse. Intermediate case c 2 — wq is the 
well studied simplest dark energy model with w = const. 
The dynamics of expansion of the multicomponent Uni- 
verse in the past depends slightly on the value of c a (see 
Fig. but depends strongly on the values of tide and wq 
related to the current epoch. This gives the possibility to 
constrain them by observational data. In this paper we 
have concentrated on determination of all parameters of 
scalar field models of dark energy. For this purpose we 
have analyzed the evolution of density and velocity per- 
turbations of scalar fields (equations (fl9"|) and (|2"0]0 and 
their influence on the evolution of matter ones. The en- 
tropy perturbations are inherent for scalar fields, so the 
dark energy pressure perturbations are connected with 
the density ones via the effective sound speed c 2 instead 
of the adiabatic one c 2 . For the classical and tachyonic 
Lagrangians it equals 1 and — w correspondingly. The 
results, presented in Fig. [5] and |H1 show that the fields 
with different c 2 and c 2 are distinguishable in principle 
if the rest of cosmological parameters is known. Unfor- 
tunately, it is not the case and we have to determine the 
parameters of scalar field models of dark energy fide, Wq, 
c 2 jointly with other cosmological parameters, the mini- 
mal set of which is: fif,, Q c dm, n s , A s , H , t, A S z- We 
restricted ourselves to the spatially flat models - that 
reduces the number of free model parameters to 9. 

For the calculations of evolution of perturbations in 
all components, the power spectra of matter density per- 
turbations and cosmic microwave background anisot ropy 
we have used the publicly available code C AMB [27], |28| , 
modified to include the presented here expressions (fT4)) - 
(Unj) and (Tll)j) - (|2"2"|) . To determine the parameters of scalar 
field models of dark energy jointly with rest of cosmolog- 
ical ones we have performed the MCMC analysis for set 
of current observational data, which include the power 
spectra from WMAP7 [4] and SDSS DR7 0, the light 
curves of SNIa [ll[, the Hubble constant measurements 
Q and Big Bang nucleosynthesis (BBN) prior (26|. It 
has been obtained that the best fitting parameters of 
scalar field models of dark energy have next values and 
lcr confidence limits: Q de = 0.72±j^, w = -0.93±£j$, 
c 2 = —0.97^0 03 hi the case of Klein-Gordon Lagrangian 
and Q de = 0.72ig;{& = -0.97tg:& ^ = -0.99±g;g? 
in the case of Dirac-Born-Infcld one. The first two pa- 
rameters, corresponding to the current epoch, are deter- 
mined well, while the last one much worse (Fig. [8]). So 
the used observational data prefer the scalar fields with 
w' > 0, which lose their repulsive properties, and predict 
the collapse of spatially flat Universe. But other models 
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with w 1 < 0, (? a > wq and even <? a = 0, are in the ler range 
of the best fitting one (Table ITTj). Therefore, the values 
of EoS parameters of the scalar fields at the early epoch 
are not determined surely enough to say which field dom- 
inates now, with receding (w' > 0) or raising (w' < 0) 
repulsion, and, correspondingly, which is the future fate 
of our Universe, eternal accelerated expansion or decel- 
erated one and collapse. We also conclude that the used 
observational data give no possibility to distinguish two 
different scalar fields (Fig. ® and Table [TTJ) , classical and 
tachyonic, but we hope that the operating and planed ob- 
servational programs will give the possibility to constrain 
better the parameters of scalar field models of dark en- 
ergy and the number of admissible evolution trajectories 
of our Universe. 
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